Horb=19100;
T=60;
Rz=6400;
Rsv=Horb+Rz;
fprintf('Rsv= %f km \n', Rsv);
t_orb=0:60:21600;
Rsv_x=Rsv*cos(2*pi*t_orb / 43200);
Rsv_y=Rsv*sin(2*pi*t_orb / 43200);
figure(1);
plot(Rsv_x,Rsv_y);
rv_x=Rsv_x;
rv_y=Rsv_y-Rz;
figure(2);
plot(rv_x, [rv_y;Rsv_y]);
V_x=diff(Rsv_x)/T;
V_y=diff(Rsv_y)/T;
figure(3);
V_sv=sqrt((V_x).^2+(V_y).^2);
plot(t_orb(2:end),V_sv);
modul_rv=sqrt(rv_x.^2+rv_y.^2);

for i=1 : length(V_sv)
    V_viz(i)=(V_x(i)*rv_x(i)+V_y(i)*rv_y(i))/modul_rv(i);
end   

figure(4)
plot(t_orb(2:end), V_viz);

c=3e+5;
nes_chast=1.6e+9;
f_dop = nes_chast*(1+V_viz/c)
k1=0;
k2=0;
for k=2 : length (t_orb)
    if Rsv_y(i)>0
        k1=k1+1;
        dop(k1)=f_dop(k-1);
        x1(k1)=rv_x(k);
        y1(k1)=rv_y(k);
    else
        k2=k2+1;
        x2(k2)=rv_x(k);
        y2(k2)=rv_y(k);
        
    end
end
 figure(5)
hist(f_dop);